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FLUCTUATION RESULTS FOR HASTINGS-LEVITOV PLANAR GROWTH 


VITTORIA SILVESTRI 


Abstract. We study the fluctuations of the outer domain of Hastings-Levitov clusters in the small 
particle limit. These are shown to be given by a continuous Gaussian process T taking values 
in the space of holomorphic functions on {\z\ > 1}, of which we provide an explicit construction. 
The boundary values W of F are shown to perform an Ornstein-Uhlenbeck process on the space of 
distributions on the unit circle T, which can be described as the solution to the stochastic fractional 
heat equation 

= —(—A) 1/,2 W(t, $) + V2 £(t,0), 

where A denotes the Laplace operator acting on the spatial component, and £(t, $) is a space- 
time white noise. As a consequence we find that, when the cluster is left to grow indefinitely, the 
boundary process W converges to a log-correlated Fractional Gaussian Field, which can be realised 
as (—A) -1 / 4 W, for W complex White Noise on T. 


1. Introduction 

In 1998 the physicists M. Hastings and L. Levitov introduced a one-parameter family of contin¬ 
uum models for growing clusters (A n )n.>o on the plane PI, which can be considered as an off-lattice 
version of discrete planar aggregation models such as the Eden model or Diffusion Limited Aggre¬ 
gation (DLA). In this paper we focus on the simplest of these models, so called HL(0), which has 
proven to be already very rich from a mathematical point of view, and has received much attention 
in recent years mmmm- 

Let D denote the open unit disc in the complex plane, and set Kq = D. At each step, a new 
particle P n attaches to the cluster K n -i according to the following growth mechanism. Fix P C C\D 
to be a (non-empty) connected compact set having 1 as a limit point, and such that the complement 
of K = O U P in C U {oo} is simply connected. Let Dq = (C U {oo}) \ Kq and D = (C U {oo}) \ K. 
Then there exist a unique conformal isomorphism F : Dq -a D and a unique constant c E M+ such 
that F(z) = e c z + 0( 1) as \z\ —> oo. We think of F as attaching the particle P to the closed unit disc 
D at 1. The constant c is called logarithmic capacity of K, and can be interpreted as the expected 
value of log \Bt\, for B planar Brownian Motion started at oo, stopped at the first hitting time T 
of K (cf. Proposition [TJ . 

Set G = F _1 , and observe that G(z) = e~ c z + 0(1) as \z\ —> oo. Let (0 n )n>i be a sequence of 
i.i.d. random variables with Q n ~Uniform[— ir, n), and set 

F n (z) := e ie "F(e~ ie "z) , G n (z) = F~\z). 

Then the map F n attaches the particle P to the unit disc at the random point e*® n . Define 
& n (z) := F\ o ■ ■ ■ o F n (z ), D n = <I> n (Do) and K n = (C U {oo}) \ D n . We say that the conformal 
map grows an HL(0) cluster up to the n-th particle, while T n = maps it out. Note that, 
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by conformal invariance, choosing the attachment angles to be uniformly distributed corresponds 
to choosing the attachment point of the n-th particle according to the harmonic measure of the 
boundary of the cluster K n -1 seen from infinity. 

With this notation, then, one has 

Pn+i — o F n j r \ o r n (ZH n ) . 

This suggest the following interpretation for the attachment mechanism: given the cluster K n . first 
map it out via T n , then attach a new copy of the particle P to the unit circle at a uniformly chosen 
point e*®”, and finally grow back the cluster K n . It is then clear that, although we are attaching 
identical copies of the particle P at each step, the particle shape gets distorted each time by the 
application of the conformal map <!>„, as shown in the figure below. 



In [25] J. Norris and A. Turner showed that, under mild assumption on the particle P, in the 
limit as n —> oo and nc —> t £ M+ the shape of HL(0) clusters is given by a disc of radius e t centred 
at the origin. Moreover, with high probability each point outside the unit disc is moved radially 
by the cluster growth. This result can be rephrased as follows: given any z G Do, <& n (z) ~ e f z as 
n —>• oo. 



Figure 1. Simulation of a HL(0) cluster with 100, 5000 and 100000 particles. 


In this note we investigate fluctuations of the map <b n around this deterministic limit. Our results 
can be divided into local and global fluctuations. 

Local fluctuations. Suppose we fix z £ C \ D, say z = e a+ia for some a > 0, a E [—7r,7r), and 
look at the limiting fluctuations of log<I> n (z:) around its mean as the HL cluster grows (n oo) 
and the point z approaches the unit disc radially (a —> 0). We prove that, provided a —> 0 slowly 
enough, these limiting fluctuations are Gaussian. Moreover, approaching the unit disc radially from 
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different angles results in asymptotically independent fluctuations. Our main result for this regime 
is the following. 

Theorem (Local fluctuations). Pick any t > 0, and let z = e ia+a for some a E [—7r, it), a > 0. 
Then as n -A oo, nc -A t and a —> 0 slowly enough, 



Moreover, the correlation between fluctuations at two different points, say z = e m+a and w = e ib+a , 
vanishes in the limit, unless the angle a — b converges to 0 fast enough with a. 


A precise statement of the above result is provided in Theorem 10 of Appendix [A] Note that the 


only Gaussian random field on T = 3D having the correlation structure defined above is the one 
given by an uncountable collection of i.i.d. A/"(0,1) random variables indexed by points on the unit 


circle. Although almost surely finite at every point, this random field is very wild, in the sense that, 
apart from not being continuous a.s., it is not even separable (i.e. it cannot be recovered by only 
looking at a countable collection of points). 

Global fluctuations. At the price of keeping z away from the unit disc while the cluster grows, 
we see that the fluctuations of logd> n become rather well behaved. In order to emphasize the 
dependence on t in this regime, and to ultimately view the limit as a stochastic process, we now 
assume nc -A 1 as n -A oo, and study the asymptotic behaviour of the map ( I> y nt \ for t E [0,oo). 

Then, by the same techniques that allow us to prove the local fluctuations result, we can show 
(cf. Theorem [lj that, for any fixed z E C \ D, these fluctuations are again centred Gaussian, with 
variance now depending on \z\ and t. Moreover, the correlation structure is sufficiently well behaved 
to enable us to prove a functional central limit theorem for log i nt i when restricted to any circle 
of the form rT := {z : \z\ = r}, r > 1 (cf. Theorem [4]). Finally, we push our analysis forward to 
obtain limiting fluctuations of log $\ nt \ viewed as a cadlag stochastic process taking values in the 
space of holomorphic functions on C \ D. Our main result is the following. 

Theorem (Global fluctuations). Let JiT denote the space of holomorphic functions on {|z| > 1}, 

and for n > 1 set P n (t, z) = ^ (log ^ ^ nt J -|_ n ^J c )- Then there exists a continuous stochastic 

process T = [Pit,-))t>o taking values in Jt? such that T n —»• T in distribution as n oo, with 
respect to the Skorokhod topology on the space of cadlag functions from [0, oo) to JP. Moreover, 
T can be obtained as the holomorphic extension of its boundary values on {|z| = 1} to the outer 
unit disc {| 2 :| > 1}. These boundary values are given by a distribution-valued stochastic process 
W = (W(t, -))t>(u formally defined in Fourier space by 



for (Afc)^, (-Bfc)fc independent collections of i.i.d. Ornstein-Uhlenbeck processes on M, solution to 



where (/?&)&, (fi'iflk independent collections of i.i.d. Brownian Motions on M. 









4 


V. SILVESTRI 


A precise statement of this result is given in Theorem [5j This provides an explicit construction of 
the limiting Gaussian holomorphic field T and, perhaps more interestingly, of its boundary values. 
Note that, since Ak(t), Bk(t) —> m-) in law as t —> oo, we have 

for ( 0 ^)^, ( bk)k independent collections of i.i.d. standard Gaussian random variables on M. This is 
(the complex version of) a well-known Fractional Gaussian Field (FGF) on the unit circle T, and 
it can be realised as Woo = (—A) - 1 / 4 kF, where A denotes the Laplace operator, and W is white 
noise on T. We remark that, in general, the FGF on T given by (— A) S W defines a true function 
only for s < —1/4, and otherwise it takes values in the space of distributions on T. It can be shown 
in greater generality (see the discussion in [ 6 ] for FGFs on M rf , d > 1) that in correspondence of 
the critical value of the parameter s , which depends on the dimension of the underlying space, one 
obtains a log-correlated Gaussian field. 


Overview of related work. Fractal patterns are ubiquitous in nature, and many attempts have 
been made to obtain a rigorous mathematical description of their formation. 

In 1961 M. Eden [9j introduced a lattice-based model of growing clusters, now called Eden model. 
The growth mechanism is as follows: start with a single site, and at each step choose one site on 
the outer boundary of the current cluster (i.e. the set of sites outside of the cluster, adjacent to at 
least one of the cluster’s sites) uniformly at random, and add it to the cluster. 

A similar growth model, so called Diffusion Limited Aggregation (DLA), was introduced by T.A. 
Witten and L. Sander [29] in 1981. In this model, at each step a new site is sampled according to 
the harmonic measure of the boundary of the current cluster from oo, and then added to the cluster. 
When, instead, one samples the new site according to the 77 -th power of the harmonic measure of 
the cluster boundary from 00 , one obtains the 77 -Dielectric Breakdown Model (DBM^), which was 
introduced by L. Niemeyer, L. Pietronero and H.J. Wiesmann [23] in 1984. Thus, the family of 
DBM, ? models interpolates from the Eden model (77 = 0) to DLA (77 = 1). 

Although these models are simple to define, they have proven very difficult to study rigorously. 
Moreover, some features of large clusters, such as the Hausdorff dimension, have been found [2, 22] 
to depend on the underlying lattice structure, which is somehow unsatisfactory. 

To overcome these problems, in 1988 M. Hastings and L. Levitov m proposed a one parameter 
family of off-lattice models, so called HL(a) models for a £ [0,oo). The case a = 0 has already 
been described in detail. When a > 0 the growth mechanism is the same, except that at each step 
the size of the new particle P n is renormalised so that its logarithmic capacity is given by 

c 


Cn — 


l$n-l( e<0 ")| £ 


( 1 ) 


If a = 2 this results in growing clusters of particles roughly of the same size, and in general ([!]) has 
the effect of attenuating the natural distortion of the a = 0 model. In m Hastings and Levitov 
argue by comparing local growth rates that the choice a = 1 should correspond to the Eden model, 
a = 2 to DLA and in general a € (1, 2) to DBM^_i. 

Although a > 0 is needed for these models to be realistic, the re-normalization 0 creates long 
range dependences which make them very difficult to analyse. In fact, to the best of our knowledge 
there are no rigorous results on (the non-regularised version of) HL(cc) models for a > 0. A first 
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regularised version of HL(a) appears in |27], in which S. Rohde and M. Zinsmeister obtained bounds 
for the Hausdorff dimension of suitably regularised clusters for a. £ [0, 2]. More recently, a different 
type of regularization was considered by F. Viklund, A. Sola and A. Turner in US!, where they 
showed that the limiting shape of regularised clusters is given by a disc for any a > 0, provided 
that the regularization is strong enough. 

The case a = 0 does not feature such long range dependences, and is much better understood. 
In {27] Rohde and Zinsmeister obtained a scaling limit for HL(0) clusters as the particle size is kept 
fixed while n —> oo. Moreover, they showed that the boundary of these limiting clusters is almost 
surely one-dimensional. More recently, in [27)1 23] Norris and Turner obtained a detailed description 
of HL(0) clusters in the small particle limit. More precisely, they proved that as n —> oo, c -A 0 
and nc —> t £ (0,oo), these clusters almost surely fill a disc with radius e l centred at the origin. 
Moreover, they showed that the ancestral tree structure within the cluster converges, always in the 
small particle limit, to the Brownian Web HU, thus providing an interesting connection between 
these two a priori unrelated models. 


Organization of the paper. We present a detailed proof of the global fluctuations result. The¬ 
orem [To] on local fluctuations is obtained by the same arguments, and we leave its discussion for 
Appendix [A] The paper is organised as follows. In Section [2] we collect some preliminary estimates 
for the basic conformal maps F. G. We then introduce in Section [3] our main tools, namely two 
sequences of backwards martingale difference arrays (10), and prove Theorem [l] on pointwise fluctu¬ 
ations. This result is then generalised in Sections [4]j5j to obtain a functional central limit theorem 
for the process (t, z) H > log <f> \ nt \ ( z ), viewed as a stochastic process taking values in the space of 
holomorphic functions on {\z\ > 1} (cf. Theorems [4] and [5]). Finally, in Section [6] we present an 
explicit construction of the boundary values of the limiting fluctuation process. These are shown 
to be given by a distribution-valued continuous process, which is rigorously defined as a Ornstein- 
Uhlenbeck dynamics in a suitable infinite-dimensional Hilbert space. We conclude the paper by 
collecting some open questions in Section [7] 


Acknowledgements. I am extremely grateful to James Norris for his guidance, support, and help 
with many of the arguments in the paper. I am also very thankful to Alan Sola for several interesting 
discussions, and to Henry Jackson for providing the simulations in Figure [TJ 


2. Preliminary estimates 

Fix a particle P as in the introduction, and, if D denotes the open unit disc, let Kq = ID), 
K = BU P. Moreover, set Dq = (C U {oo}) \ Kq, D = (C U {oo}) \ K. R follows from the 
Riemann Mapping Theorem that there exists a unique conformal maj0 F : Do —> D such that 
F(z) = e c z + 0(1) as \z\ —> oo, for some constant c £ M. We assume that the particle P is regular 
enough so that F extends continuously to the boundary of Dq. Set G = F ~ 1 . 

Conformal maps are very rigid, and simply from the definition of F. G we can deduce the following 
properties: 

(PI) |T(^)| > \z\ for all 2 £ Do, |G(^)| < \z\ for all z £ D, 

^Throughout this paper, by conformal map we mean a conformal isomorphism. 
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(P2) there exists a constant C > 0 such that | F(z)\ < C\z\ for all z G Do and \G(z)\ > \z\jC for 
all z G D. 

Indeed, (PI) follows trivially from Schwarz lemma, while (P2) is a consequence of the prescribed 
behaviour at infinity for F. G. 


Notation. Throughout the paper, C denotes a finite, positive constant which can change from line 
to line, and which is either absolute or only depends on the maps F. G. Whenever this constant 
depends on other parameters, say a, (5 ..., we make it explicit in the notation by using C(a, (3 ...). 

In order to obtain finer distortion estimates for the maps F. G it is often useful to relate the 
logarithmic capacity c to geometric properties of the particle P. Assume the following: 


Assumption 1 . There exists 5 > 0 such that 

P C {z G C : \z- 1| < 5}, 1 + 5 eP, P = {z:zeP}. (2) 


We regard 5 as measuring the diameter of the particle P. 

Remark. Assumption [l] is in force throughout the paper. 

The following result appears in [25] (cf. Proposition 4.1 and Corollary 4.2 therein). 


Proposition 1. There exists an absolute constant C > 0 such that, for all z G D: \z — 1| > 25, it 
holds: 


log 


G(z) 


+ c 


< 


Cc 

\z — 1 | ’ 


d 1 

dI log 


G(z) 


< 


Cc 


- II 2 ' 


Moreover, 

for 5 small enough. 


5 2 3<5 2 

— < c < - 

6 “ “ 4 


( 3 ) 

( 4 ) 


In light of Q we use c and 5 interchangeably, and all statements are intended to hold for c, 5 
small enough. Combining ©> and Q we deduce the following improved bound. 


Corollary 1. 

holds: 


There exists an absolute constant C > 0 such that, for all z G D: \z — 1| > 25, it 



+ c 


z + 1 


z - 1 


< 


Cc 3 G\ z \ 

I z — 1| 2 ' 


( 5 ) 


Proof. We follow the proof of [25] . Proposition 4.1. Let u,v denote the real and imaginary part 
of log — z 1 respectively, so that they are harmonic functions on D. Then by optional stopping 
u(z) = — E(log |-Bt|) < 0, T being the first hitting time of K for a Brownian Motion B starting 
from 2 . Introduce the particle Pi D P defined by P\ = {z G Dq : \ z z f\ \ < r}, for r = <5/(2 — 5), 
and set D 1 = (C U {oo}) \ (DU P\). Then the unique conformal map G\ : D 1 —> Dq satisfying 
Gi(oo) = oo and G\ (oc) > 0 is given by 


Gi{z) 


z{+iz — 1 ) 
z — 7 


for 7 


1 — r 2 
1 + r 2 


Set F\ = G l 1 , and A = {z G dP\ : \z\ > 1}. Then G\(A) = {e *’ 9 : |$| < ■do} with do = cos 1 7 . 
Moreover, u o F\ is harmonic and bounded on Dq and, using that ^ 0 l 7 i)(e* , 9 )d?? = —c, 
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the optional stopping theorem yields 


{uoF^z) = -c+ — 


(u o FiXe^) Re 




2 e 


I'd 


z — e 




It follows that 

(u o Fi)(z) + c + Re 


2 c 

z- 1 


= £/ («"?!)(«*') 


2tt 




Re 


2 e 




z — e 




di9. 


— Re 


z - 1 /J 


dtf. 


Using 0 , then, we find 1 — cos??o = 1 — 7 x S 2 /2 < Cy/c, from which 

G|z|-y/c 


Re 


2e 




z — e 

I'd 


I'd 


— Re 


2 

z- 1 


2 |z||l — e il?0 | 

~ \z — e 1 ® 11 z — 11 — 


z — e l,J | |z — 1 | 


< 


Since (u o .Fi)(e”j < 0 for all i? in the integration range, we find 


(u o F 1 )(z) + c + Re (y—~ y) 


< 


C\z\yfc 


dist(z, Gi(A)) 2 


± f (uo Fi)(e^)M 

J\0 1<1? 0 


C\z\y/c 

dist(z, G\(A)) 2 ’ 
C|z|c 3 / 2 


Now take Gi(z) in place of z, to get 
2 c 


u(z) + c + Re( j^i) < “P)+ c + r XgI(7Wt) 

C|z|c 3 / 2 


+ 


Re 


2 c 


Gi(z) - 1 


dist(z, Gi(^4 )) 2 

(^t) 


— Re 


+ 


2 c|Gi(z) — z | 


( 6 ) 


“ dist(Gi(z),Gi(^4 )) 2 |Gi(z) - l||z- 1| ’ 

where in the second inequality we have used that |G(z)| < \z\. Using the explicit expression for Gi, 
one shows that |Gi(z) — 1| > |z — l |/2 and |Gi(z) — z| < Cyfc\z\ for 5 small enough, from which 


2c|Gi(z) — z | 


< 


C|z|c 3 / 2 


(7) 


|Gi(z) — 111z — 11 ' |z — 1_ 

Moreover, reasoning as in [25] one shows that there exists an absolute constant C\ such that 
dist(Gi(z), G\(A)) > |z — l|/Gi. Putting this together with ([h]) and Q we finally obtain 


u(z) + c + Re ^ i ) 


< 


C|z|c 3 / 2 
|z — 1| 2 


for all z E D such that |z — 1| > 25. 

Now note that u(z) + c + Re (zzri) and v(z) + Im (^-) are the real and imaginary part of the 
holomorphic function z i —> log + c + on D n {z : |z — 1| > 2<5}. It then follows by Cauchy’s 
integral formula that 

C|z|c 3 / 2 


vfr>(z) +Im( 2 ^ i )^ = V ^u(z) + c + Re ( 7 ^-y)^ 
Finally, using that | v(z) + Im (^rj) | —> 0 as |z| —> 00 , we get 

2 c 


< 


2 c 


v(z) + Im y- - -j < 

which concludes the proof. 


r v( v(z + s(z — 1)) + Im ^ 


z + s(z — 1 ) — 1 


- II 3 


|z — l|ds < 


C\z\6 


3/2 


|z- 1| 2 ’ 

□ 
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We combine Proposition [T] and Corollary [l] to obtain corresponding estimates for the function F. 
which are collected below. 

Corollary 2. There exists a constant C > 0 such that, for all z G Dq with \F(z ) —1| > 26, it holds: 


| F(z) - e c z | < 


Cc\z\ 
\z\ - 1 


log 


F(z) 


z + 1 


— c 


z- 1 


< 


C-c^/ 2 \z \ 2 

(M-i) 3 ' 


Proof. For the first inequality, note that ([3]) readily implies that 

Cc|z| 


| G(z) — e c z | < 


\z - II 


( 8 ) 


for all z E D : \z —1| > <5, and 6 small enough. Therefore, by setting w = F(z) and using (P1)-(P2), 
we obtain 

Cc\w\ , C'c\z 


| F(z) — e c z | = e c | G(w) — e c w\ < 


< 


(9) 


\w — 1| — \z\ — 1 ’ 

for all z E Dq such that | F(z) — 1| >26 and <5 small enough, as claimed. 

For the second inequality, note that, since \F(z) — 1| > 26 by assumption, we can use ([5]) to get 

1 1 


. F(z) z + 1 
log-c 


2 - 1 


< 


, G(w) w + 1 
log + c 


w 


w — 1 


+ 2c 


w — 1 G(w) — 1 


^ Cc 3 / 2 \w\ + 2c\G(w) — w\ < Cc 3 / 2 |z| ^ 2c\G(w) — w\ 


\w — 1| 2 \w — l||G(u;) — 1| (|z| — l) 2 (\z\ — l ) 2 

Moreover, it follows from Q that 

| G(w) - w\ < | F(z) - e c z | + (1 - e _c )|u;| < + Cc\z\ < 


for c small enough. Putting all together, we end up with 

C(?! 2 \z 


, F(z) z+1 
log- c 


z z — 1 

for c small enough, as claimed. 


< 


(M-i ) 2 


x i y/c\z\ \ < 2Cc 3 / 2 |z| 2 


z - 1 


(M-i ) 3 


□ 


3. POINTWISE FLUCTUATIONS 

In this section we prove that, fixed t > 0 and 2 G C\B, the fluctuations of log + i nf i (z) around its 
mean are given by a complex Gaussian random variable, whose variance is independent of Arg(z). 

Notation. Throughout the paper A f(/a,a 2 ) denotes the Gaussian distribution on M with mean /x 
and variance a 2 . 

Our main result is the following. 

Theorem 1. Fix any t > 0. Pick a € [—7r, 7t) , a > 0, and let z = e m+a . Define v 2 (o) := 
log 1 ^ 1 e _ e l 2 a ’ . and let F a (t,e ia ) be a complex Gaussian random variable with i.i.d. real and imagi¬ 
nary part, distributed according to Af(0, v 2 (a)). Then it holds: 

>g * L ^ _ [nt\cj —> F°{t,e ia ) 


in distribution as n -+ oo, c —> 0 and nc —> 1. 
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We prove Theorem [l] for a < 1, which we assume without further notice. This entails no loss of 
generality (see discussion in Section[5]), and it has the advantage of slightly simplifying the notation. 

Remark. From this point onwards, with the exception of Appendix [A] whenever we write n —> oo 
or c -) 0 we mean that n -A oo, c —> 0 and nc —> 1. 


Our main tool for the proof of Theorem [l] consists of two sequences of backwards martingale 
difference arrays, that we now define. Note that Ln 1 — > 1 for all z £ D o, so log ——— defines a 


holomorphic function on Dq. We fix the branch of the logarithm by requiring that log ^ ut J 
as z —> oo. For a > 0 and a £ [—7r,7r) as in Theorem [lj then, we find: 

0 ia+a\ l nt l i „ Z7. „ „ jp, . („ia+i7\ L n fl FVo - *0 


(*) 


. $lnti(e ia+a ) F k oF k+1 o 

lQ g -- = Z. lQ g 

k =1 


fitlO'-oF|„j(e“«) A. ° g 


where (a) := F k+1 o • • • o F ^ nt j (e m+<T ). Moreover, if ,F k . n := cr(Q k , 0 fc+1 ... 0 n ), then 


E ( log ■ 




e iBkZ hnt}( a ) 


^k+l,[nt\ ) = ^ / lo S 


F(e 


—id > 


We therefore set 

1 


x: 


k, [nt\ 


(«) = 




log 


ne- l ^Zl [nti (a)) 

e ~ iekZ k,[nt \( a ) 


— c 


^-J (a)) »_ , im ,„ g ^> 


e |z|->oo £ 


= c. 


k , |_ra£j 


1 






( 10 ) 


where Arg(z) £ [—7r, 7r). The above computation then shows that (a))j<[ n tj and (XZ^j (a))fc< 

form sequences of backwards martingale arrays with respect to the same filtration (J^ k,n)k<\nt\ ( see 
[3], Section 35 for the definition of martingale arrays). Moreover, 


\nt\ 

£ 

k =1 


1 / <I>i t \(e ia+CT ) \ 

{ X k,[nt}( a ) + iY k,[nt\( a )) = ^ lo g [n e ia+a - \- Ut \ C ) ‘ 


Remark. Throughout the paper we identify C with M 2 , and often refer to complex random variables 
as random vectors and vice versa, depending on which point of view we seek to emphasize. 


The following result appears in ED, Corollary 2.8, as a Central Limit Theorem for (forward) 
martingale difference arrays (see also [3], Theorem 35.12). It is straightforward to adapt the proof 
to backwards martingale difference arrays, so that we have the following. 


Theorem 2. Let ( X k n )i<k<n be a backwards martingale difference array with respect to ' k n = 
&[ x k,n ■ ■ • n )■ Let S k , n — ^ b —b Assume that. 

n 

(I) for all rj > 0, X k n \{\X k ^ n \ > rf) — > 0 in probability as n —> oo, 

k =1 
n 

do £ X k n -A s 2 in probability as n -A oo, for some s 2 > 0. 
k=1 

Then S rhn converges in distribution to Af(0,s 2 ). 

Note that the above theorem is concerned with scalar random variables, while we have 2- 
dimensional vectors. In order to reduce to the scalar case, recall that by the Cramer-Wold 
Theorem (cf. 0, Theorem 3.9.5) it suffices to prove convergence in distribution of all linear 
combinations of the vector entries. To this end, pick any a, (3 £ M, and note that by linearity 
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(aX^ nt j(a) + \nt\^ a ))k<\nt\ a S & i n a backwards martingale difference array with respect to 

the filtration {^k,\nt\)k<\nt\- We are going to apply Theorem [ 2 ] to this linear combination. To this 
end, we collect here some estimates for (Xf (a)) and (Y£ y n t\( a ))- Since a and o are fixed, we 
omit them from the notation throughout this section. 


Lemma 1. There exists a constant C > 0 such that, for c small enough, it holds \X k \ nt \ \ < C/yfc, 
\Yk,[nt\ I < C/^c for all n>l,k< [nt\. 


Proof. It follows from (PI) that log | | > 0, from which \X k | > —yfc. Moreover, 

(P2) gives \X k \ nt \\ < Cfy/c-V y/c < 2C / y/c for some constant C > 0 and c small enough. Finally, 
since Arg( F ^ iefc ^’ [ri | j) ) E [—7r, t r), we have \Y k ^ nti | < tt/\/c- □ 


J k , |_ntj 


Much better estimates can be obtained by restricting to a certain good event, which is shown in 
to have high probability for large n. The following result identifies this event. 


Theorem 3 (|25j. Proposition 5.1). Fix a positive integer m and a constant £ > 0. For n < m 
define the events 

E n {e) : = {|e _cn <l> n (z) — z\ < ee 6e for all z : \z\ > e 5e } 

n {\e cn T n (z) -z\< ee 5£+cn for all z : \z\ > e cn+4e } , 
and set E(m,e) := H™=i E n (e). Then it holds 

P (E(m,e) c ) < C(m + £~ 2 )e~ £3 ^ Cc ^ 

for some constant C > 0. In particular, by setting m = |_5 _6 J and £ = <5 2 / 3 log(l/<5), one obtains 
that 5~ k P(E(m,£) c ) —> 0 as 5 —> 0, for any k > 0. 


We refer to E{m, e) as the good event. 

Remark. Without further notice, we take m = [<5~ 6 J and e = <5 2,/3 log(l/<5) as in the last part of 
the above theorem, so that e><5 and n x 5~ 2 <C m. 


Lemma 2. Assume that a 2> <5. Then there exists a constant C{t) > 0 such that, for n large 
enough, on the good event E(m, e) it holds 


mj Sj {|jr i , LnlJ |v |y fc ,L„. j l}<c ( t ) A 


(ii) 


Proof. For any k < \ nt\ we have 

max {| X k 1 n t\ |, \Y k j nt\ I} — y= 

V c 

1 



log F(e- i0 *Z kAnti ) - log(e- i0 ^ fc)LntJ )| + c) 

sup rrr) • \ F ( e ~ l&k z k,lntj) - e~ ie>k Z k ^ nt j | + c 
|£|>e CT Isl ' 


where the last inequality follows from the mean values theorem, and the fact that \F(e i<s>k Z k int|)| > 
\e~ lQk Z kt y nt ^ | > almost surely by (PI). Now note that \F(e~ l&k Z k ^ nti ) - 1| > \Z k ^\ - 1 > 
o ^$> 25, so by Corollary [2] we have 
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for z = e l&k Z k \ nt \. Moreover, since we are on E(m,e ), there exists a constant C(t) depending 
only on t such that e a < \Z k \ nt \ < C(t). This, together with (12), yields 


max < \X, 


k< \nt\ 

as claimed. 




V \Yi 


k,\_nt] I [ — 


< 


1 /'2Cc\Z k ^ nt ^\ 


\/~C V l^fc,Lrctjl 1 


+ C < 


C(t)y/c 


V~c 


e a — l +^^ 2C W- — 


□ 


We now make use of the above bounds to prove that the backwards martingale difference array 
(^fc, \nt\ )fe< \nt\ i with x k,\nt\ ’■= oiX k ^ nt j + PY k ^ nt J, satisfies Assumptions (I)-(II) of Theorem [2j In 
doing so, we provide an explicit formula for the limiting variance. 


|ntj 


Lemma 3. Assume that a 5. Then for all r/ > 0 it holds E ^Lntj^d^fc.Lntjl > rf) 0 in 

k =1 

probability as n —> oo. 

Proof. For any e > 0 we have: 

> e ) ->”) - HkS 14 ’ 1 "' 1 ) 

= (H I; £ ) C ) + ^ ^ Kl«J I; E (m,s)j . 

The fact that the first term in the r.h.s. converges to zero as n —> oo follows from Lemma [l] and 
Theorem |3j while convergence to zero of the second term is a straightforward consequence of Lemma 

E □ 

We now concentrate on Assumption (II). The first step consists in replacing condition (II) with 
a more convenient one, involving conditional second moments. The following result shows that, 
provided a is large enough with respect to c, this is allowed. 

|ntj 

Lemma 4. Assume a 3> yfb, and that ^^E(A^ . t , \^ k +i ,[nt] ) s 2 in probability as n —> oo for 


k= i 


[nt\ 


some s 2 > 0. Then also |ntj s2 probability as n —>• oo. 


k =1 


Proof. Let M k>LntJ := X%^ nti - I^+Untj)• Tt is readily checked that (M fcj |_ nt j) fe <|_ n tj is a 

backwards martingale difference array with respect to the filtration {&k,\nt\)k<\nt\ • We aim to show 

|ntj 


that for any rj > 0 it holds P (IE M k , V nt\ 


k =1 


> T] 


\nt\ , j- [ntj 

5Z M ^tj| >t?) < -jE( 

fc=l ' VL fc=l 


r L nt J n 2 

k, nt\ 


0 as n -» oo. Indeed, 


LntJ (+) x LntJ 


iE E K t »<j) < 4 E e «m) 


k =1 


fc=l 


LntJ L n d 

= ? E ; £ (">. >T) + -3 E ■ 

' fc=i ' fc=i 

Above, (*) follows from the general inequality E((X — E(X)) 2 ) < E(A 2 ), which we apply to each 
term with respect to E( • \.Y k i nt i ). The fact that both terms in the r.h.s. converge to zero as n —> oo 
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is now a consequence of the bounds for \X k \ nt \ and \Y k and hence for l-'t'*. i nt j |, obtained in 
Lemmas HH2 □ 

In light of the above result, it remains to compute the limit in probability of 
\nt\ [ntj |n*J 

y ^( X k,\nt\ \^k+l,\nt\ ) = « 2 y : |ntj l^fe+l,LntJ ) + Z^ 2 ^ E (*fc 2 ,[ntJ l^fc+l.lntj ) 

/c=l k =1 /c=l 

[ntj 

+ 2 a/I E ® (^fc, [ntj ^ fc, [ntj I ^fc+l, [ntj ) * 

fc=l 

and prove that it coincides with (a 2 + fi 2 )v 2 {a), where v 2 (a) is the limiting variance introduced in 
Theorem [l] The following result shows that, in fact, it suffices to compute the limit of the first term 
in the r.h.s. above. 


Proposition 2. It holds |^fc+l,LntJ ) \^k+l,[nt\ ) J an ^^^k,\nt\^k,\nt\\Z'k+l,\nt\) ~ 

0 almost surely for all k < [nt\ . 


Proof. All equalities in this proof are intended to hold almost surely. Introduce the holomorphic 
function f(z) := ^ (log — c) defined for \z\ > 1, so that 


f(e- ie *Z k , lnti ) 


1 2 


^■k, [nt\ Y k> |^ nt j + 2 iX k , ^ nt j I ^ [ntj 


Taking conditional expectations both sides, we find 


1 C n r 12 

_ / f ( r —i$ ry \ _ H7/V'2 

2t r ./ 


/( e * di? — E(A fc ^ nt j|^ fc+1 ^ rit j)—E(y fc ^ nt j|^ fc+1 ^ nt j)+2iE(X fc ^ rif jI / fc ^ rit j|^ fc+1 ^ nt j). 

(13) 

On the other hand, / being holomorphic, Cauchy’s integral formula yields 

2 


1 

27T 


/(e-‘% L ntj) 


1 2 


d$ = -—: 


J. f Zk, [nij \ 

2 d z 

— = lim 

tf ^ k , |rat| \ 

; V * )\ 

Z |z|—H) 

. V * ) \ 


= 0 . 


Gong back to (13), this implies that both real and imaginary part of the r.h.s. must vanish almost 
surely, which is what we wanted to show. □ 


Proposition [ 2 ] already shows that the limiting Gaussian vector J ra (t,e ia ) must have i.i.d. entries. 
It remains to compute the limiting variance, that is to show that for all r] > 0 it holds 


LntJ 


YMXl\nt\\^k+Unt\) ~ ' y2 ( cr ) 


k =1 


> T] 


as n —> 00 . To this end it is clearly enough to work on E(m,e), the advantage being that on this 
event we have 

for all k < [nt\, as it follows directly from the definition of E{m,e ) as long as a e. Our strategy 
is then to replace each Z k \ nt \ by its deterministic approximation, and show that, provided a is large 
enough with respect to c, this does not affect the limiting variance. 
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Recall that the Poisson kernel for the unit disc D is given by P r (i?) = Re( ^^ e w ) for r < 1, and 
that the function re™ t —> P r ('d) is harmonic in D. Moreover, given any continuous function / on 
T = cMD>, its harmonic extension Hf inside D is given by Poisson’s integral formula 

(■ Hf){re id ) = ^f MV ~ t)f(e u )dt = (P r * f)(d). 

We denote by Q r ($) the harmonic conjugate of P r (d) in D, i.e. Q r (i?) = Im() ■ 

Lemma 5. Assume a e. Then there exists a constant C(t), depending only on t, such that on 
the event E(m, e) we have 

ml^k+UntO + C-^ (P e _._ (L „ t j-fc )c (tf)) 2 dtf 

for all k < [ nt\, and n large enough. 


C(t)ce 

— (a + ([nt\ — k)c) 3 


This result follows by a more general one, namely Lemma [7] in the next section, and the proof is 
therefore omitted. Assume now that a 3> yfe. Then we deduce from Lemma [5] that on E(m,e) it 
holds: 


|ntj 


\nt\ 


\nt.\ 


'^2^( X k,[nt}\^ r k,[nt\) 2vr / [nt \c < ^2 / , (\nt\ _ m c \ 

k =1 k =1 J ~ 7C k=l V VL J 11 


< C(t)e 


fa+\nt\c ^ 


„ C (t)e( 2 


2a 2 2(a + [nt\c) 2 


as n —> oo. Note that yfe x 5 1 / 3 (apart from logarithmic corrections), so the assumption o 3> yfe 
is stronger than the previous one cr S> yfS. 

In conclusion, we have shown that that, provided a yfe, the limiting variance is given by the 
deterministic expression 


lim ( — 

n—> oo 


( l nt \ /*7T \ PCr-\-t 1 /*7T 

— ^ / (P e -.-( L n t j-/= )c (i?)) 2 di?- Infjcj = — (P e -x (id)) 2 dddx — t 

t — 7T / J G J — 7T 


p<T+t ra+t 

(P e -x * P e -x)(d)dx — t = 


1 + e~ 2x , , 1 - e" 2 ^ 2 


1 — e 


- 2 ® 


d.x — f = log ■ 


1 — e 


—2a 


= v t , 


where we have computed the inner integral by mean of Poisson’s integral formula. Finally, if a > 0 
is kept fixed as n —> oo the assumption a yfe is trivially satisfied, so this concludes the proof of 
Theorem [I] 


4. The fluctuation process on C(T) 

Having a pointwise convergence result, it is natural to ask if this can be extended to obtain 
convergence of random fields. Recall that T = {|z| = 1} denotes the unit circle, and let C(T) 
denote the space of continuous functions from T to C, equipped with the supremum norm 

IMloo = sup \x(e l ®)\ . (14) 

$(=[— 7T,7r) 

Moreover, let D[0,oo) denote the space of cadlag functions x : [0, oo) —> C'(T). The goal of this 
section is to prove the following result. 
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Theorem 4. Fix any a > 0, and let F° denote the C(T)-valued cadlag stochastic process defined 
by 


FZ(t,e ia ) = ±=( log 


$lnt\(e ia+ °) 


a icL-\-cr 


— \nt\c 


for e m E T and t > 0. Then there exists a continuous zero mean Gaussian process F a : [0,oo) 
C'(T) whose covariance structure is given by 


Cov(J^(t,e“)) =v 2 t (o) 


T 0 
0 1 


Co v(F a (t,e ia ),F a (s,e ib )) = 


c s ,t(p,a-b ) c Si t(<r,a-b) 
-c Sit (a,a-b ) c s t (a, a - b) 


where v^(a) = ct,t(&, 0), and for s < t 

/ i _ g—2 a—(t+s)+ia 

c s , t (a,a) ■■=Re^log i _ e _ 2a _ {t _ s)+ . a 


c Sit (a,a) = Imf log 


^ _ e -2o--(t+s)+ia 
^ g—2(7— (t— s)+ 2 Q! 


such that F° —> F a in distribution as n —»• oo ; m the sense of weak convergence of probability 
measures on the space D[0,oo) equipped with the Skorokhod topology. 


Note that the limiting process is rotationally invariant in the spatial coordinate, as one expects 
from the rotation invariance of the original model. The rest of this section is devoted to the proof 
of the above result. 


Fix any cr > 0. It is trivial to check that Ff belongs to Z)[0, oo) for all n > 0. Since D[0,oo) 
equipped with the Skorokhod metric is a complete separable space, it follows by Prohorov’s theorem 
that Ff — > F a weakly if and only if the finite dimensional distributions (FDDs) of F° converge to 
the ones of F a , and (F%) n >o is tight (see m, Lemma 4.3). 

4.1. Convergence of finite—dimensional distributions. The following result is a direct conse¬ 
quence of the discussion in [H]. 

Lemma 6 (|33 lJ). Assume that: 

(i) for any t > 0 the family of probability measures of (F£(t, -)) n >o on C'(T) is tight, and 

(ii) for any 0 < t\ < ... < tM and any —n < a± < ... < aju < tt, M E Z+, it holds 


( F%(ti,e iai ) 

... F% (t M , e iai ) \ 


/ F a {h,e iai ) . 

• • F°(t M ,e ia i)\ 

\F%(ti, e iaM ) 

... FZ(t M ,e ia ™)j 


yF a (ti, e*“ M ) . 

.. F°(t M F iaM )) 

distribution as 

n —> oo. 




Then the FDDs of F° converge to the ones of F a as n —> oo. More precisely, for any M E N + 
and any 0 < t\ < ... < tM, the random vector (F%(ti, ■),... ,F% (tM, ■)) * n C(T) RI converges in 
distribution to (F a {t\, •),... ,F a (tM , •))• 


We start by showing that (ii) holds. To this end, let us again reduce to the scalar case by 
considering linear combinations. Recall the definition of the random variables X% (■), Yff (■) 
given in (10). For (aij)i<ij<M in we look at the weak limit of 


\nt M \ 

^ a ljFn (ti , e mj ) = V] V] aij 


1,3 


k =1 l,j 


\nt M \ 


■^k,\nti\i a j) + ^k,[nti\( a j) < \FU\) ~ ^k,\nt M \ 


k= 1 
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where we have set 




k, \nt M \ 




X k, |nt,j (%) + iY k,\nti\ (°i) < L^J ) 


( 15 ) 


for k < [ntM\ • Then {^k\nt M \)k<\nt M \ a backwards martingale difference array with respect 
to the filtration {^k,\nt M \)k<\nt M \- We can therefore apply Theorem [2] to show weak convergence, 
provided that Assumptions (I) and (II) are satisfied. Note that all terms in the r.h.s. of ( |15[ ) satisfy 
the estimates of Lemmas [I] and [2j from which one can easily show, reasoning as in Lemma [3j that 

Assumption (I) holds. Furthermore, Lemma[4]is still in force, from which we conclude that, provided 

\nt M \ 

a VS, the limiting variance is given by the limit in probability of E E ((-*A.<«j) 2 I ^+1, [nLvf J ) ‘ 

k= 1 

We now focus on the computation of this limit. Expand the square and use linearity to see that the 
above sum equals 


\nti\/\\nt r \ r 

yy a ij a rs ^2 IE ( dj ) + iY^ (a s ) + iY^y ntr j (a s )^ 

l,j r,s k =1 L 

It follows that it suffices to compute the limit in probability of 


Y'k,\nti\ vpltj- 


|nsj 

(x*, L „ tJ (») + («>) (vj lniJ (f>) + .Y£ L „j 

fc=l 




for arbitrary a,f) 6 [—7r, 7 t) and 0 < s < f. Moreover, by rotational invariance we can set 6 = 0 
without loss of generality. The following result simplifies the computation. 


Proposition 3. Almost surely, it holds 

^( X k,[nt\ ( a ) X k,[ns\ (0)l^fc+l,|ra*j) = {a)Y£^ ns j (0)|LriiJ ) 

^ X k,\nt\ ( a ) X k,\ns\ (0)|^fe+l,LntJ ) = -E (^LnsJ ^) X k,\nt\ { a )\^k+l,\nt\ ) 

for all k < |_nsj . 

Proof. The result follows by the same arguments used in the proof of Proposition [2j considering 
now [f{er iQh Z a k Vnt j (a)) + /(e _<e *Z£ LnaJ (0))] 2 in place of [/(e" i0fc Z fc)Lnt j)] 2 . □ 

It remains to compute the limit in probability of [ntl ( a ^ X k Insl (0)|'^fc+i,|ntj) an d 

E* E W, L ntj( «)nVj(°)l^+L L ^j). P rev i° us section, we do this by approximating by a 

deterministic quantity. 


Lemma 7. Assume a S> e. XTien t/iere exists a constant C(t), depending only on t (t > s), such 
that on the event E(m,e) the following hold: 


^ X k,\nt\ ( a ) X k,\ns\ (0)l^fc+l,|n£|) + c ~ 
^( X k,\nt\ ( a ) Y k,[ns} (°)l^fc+l,|nij ) + c “ 




I?)P e _ CT _ (Lri3 j_fe) c (l?)dl? 

'd'jQ f>—& — (L ns J —k)c (i?)di? 


C(t)c£ 

~ ( <T +(L n ' s J — k)c) 3 
C(t)c£ 

~ (cr+(LnsJ — k)c) 3 


for all k < and n large enough. Above Qri'd) = Im( |+(^ ), r < 1, denotes the conjugate 

Poisson kernel. 
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We discuss the proof of Lemma [7j from which Lemma [5] also follows by setting a = 0 and s = t, in 
Appendix [bJ As a consequence of the above result we obtain that, after further assuming a » y/e, 
on the event E(m,e) it holds: 


< 


LnsJ [nsJ pn 

E(Aj^ | nt | (a)X^ | ws | (0)|#'fc +1 j wt |) + [nsjc — — / P e -v-(\nt\-k)c(a — i?)P e -a-(Lnsj-fe)c('f?)dd 

k =l 7 r / c =i , '- 7r 

L ns J __ /-cr+lnslc j __ / I 

.. f — - .7 , 1 ^0 

X 


C(t)ce 


< ^ ( CT + (L ns J -k)cf 


<-(T+\ns\c J 

< C(t)e j -3 = C(t)e 


\2ct 2 2(ct + |_nsjc) 2 


as n —> oo. This in turn implies that 

|nsj 


L ns J / L ns J pTT 

J™ Y1 E ( X k,[nt\ ( a ) X k,\ns\ (0)l«^fc+l,Lntj) = \ P e—(L»*J-*)c(a “ l?)^ e -a-(Ln S j-fc)c(l?)dl? - |nsj C 

°° fc=l °° ^ ^ fc=i •'~ 7r 


1 

2 tt 


r<T+s /"n- 


P e -(t- s )-x(a - #)_P e -x(l?)d$ — S = 


r+^T 

J(T+^ 


P e - 2 x(a)dx — s = log 


^ _ g-2cr—(t+s)+ia 


^ _ g—2 cr—(t—s)+ia 


as claimed. Similarly, 

[nsj 


l_j -j rcr-\-s rir 

}^J2 E ( X llnt]( a ) Y k,lns](°)\^k+l,lnti) = / / P e -(t-.)-*(a - 0)Q e -*(0)d0 

/ i ^ cr «/ —7r 

^ _ g—2cr— (t+s)+za 


fc=l 


/cr+ 


1 2 
t—s 


Q e - 2 X (a)dx = Arg 


^ g—2cr— (£—s)+ia / 


This concludes the proof of (ii) in Lemma [ 6 j It remains to show that (i) holds. 

Pick any t > 0. Suppose we could show that: 

(a) for all v > 0, there exists M v > 0 and no £ N+ such that sup P(||J r ^(0, ^Hoo > M„) < z/, 

n>no 

and 

(b) for all v > 0 it holds 

lirri lirn supP( sup \J 7 °(t,e m ) — P°(t,e lb )\ > v) = 0. (16) 

T ? _> ' 0 n—too \ \a—b\<r/ J 

Then it would follow from Theorem 7.5 in [Ij that the sequence of probability measures on C(T) 
associated to (P"(t. -))n>o is tight. 

Since P^(0,-) = 0, (a) is trivially satisfied. The fact that also (b) holds follows by the same 
reasoning as in Section [4.2| below, and the proof is therefore omitted. 


4.2. Tightness. To conclude the proof of Theorem [4] it remains to show that the family of proba¬ 
bility measures associated to (P%)n >o on Z?[0,oo) is tight. Corollary 7.4 in [10] provides a charac¬ 
terization for tightness in this space. In particular, it tells us that if: 

(a’) for all u > 0 and all t £ [0, oo) n Q, there exists a positive constant M = M(u. t) such that 

limsupP(||J^(t, -) ||oo > M) < v , 
n—¥ oo 

and 

(b’) for all v, T > 0 there exists a positive constant r/ = rj{y, T) such that 


limsupP(w(J r ^,r/,T) > v) < v , 

n—¥ oo 
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where uj denotes the modulus of continuity on D[ 0, oo), then the desired tightness follows. We start 
by showing (b’). It is clear that it suffices to restrict to the event E(m,e). We are set to show that 
for any u,T > 0 it holds 


limsupP sup e ia ) — (s, e lot )\ > u; E(m,e) 

n^roo y s,te[0,T],|t—s|<77 , 

aG[—7r,7r) 


0 


as r\ —> 0 . 


(17) 


It is convenient to switch to logarithmic coordinates. Following the notation introduced in 
Section 5, we set D$ = {w £ C : Re(u>) > 0}, D = {w £ C : e w £ D}, and let F be the unique 
conformal map from Dq to D such that F(w) = w+c+o( 1) as Re(u>) —> oo. Moreover, let G = F , 
so that G(w) = w — c + o(l) as Re(iu) —> oo. Finally, for all k < n set F k (w) = F(w — *©&) + iQ k , 
G k = F^ 1 and = Fi o ■ ■ ■ o F k , f k = Then F£(t, t?) = ^ (ipKl + cr) - id - ° - [nt\cj 
and, assuming t > s without loss of generality, we have: 


sup | T u n (f, e ia ) - (s, e ia ) \ > v ; E(m , e) = 

s,te[0,T],|t—s|<?7 J 

Q;G[—7T,7r) 


sup 

\ts\<T) 

z:Ke(z)=a 


®[nt}(z) -z - [nt\c] - ($ Lns j(z) - z - [nsjc 


(18) 


> vyfc .; E(m , e) 


Note that on E(m,s ) we have z = F| ni j (w) = r'i ris j('u/) for some w £ D\ n t\ and w ' £ D\ ns \ such 
that | Re(w) — cr — |_nijc| < e, | Re(w/) — a — |_ras_|c| < e. Moreover, since |t — s| < r/, it must be 
| w — w' | < [nt\ — L ns J + 4e < 4?/ as long as n is large enough. For any k > 0, if w £ D k define 
M k (w) = ffc(re) — w + kc. Finally, let 

T : = inf [k > 0 : £ ^ D k for some £ such that Re(£) = cr + k — e}, 


and note that on E(m,s) we have T > [nt\ f° r all t < T. It follows that the r.h.s. of (18) is 
bounded above by 


sup 

„ \t~s\<v 

w :| Ke(w)—a— |_ntjc|<£ 
w':\ Re(it/) — cr— [nsjc|<£ 
\w—w'\<4rj 


r |ntj (w)-w + [nt\ c) - (r [ris j (w’) - w' + |nsj c 


> v\fc ; E(m, e) 


< P sup 

\ | t — s |<?7 

w :| He(w)—a— [ntjc|<e 
w':\ Re(w')—cr— [nsj c\<e 
\w—w'\<4r) 


M lnt\ATi w ) ~ M [ns ja r(V) 


> Ua/c 


We control the above probability by mean of a 3-dimensional version of Kolmogorov’s continuity 
theorem (cf. (7], Theorem 1.6). To this end, we show the following. 

Lemma 8. Fix any s, t £ [0, T] with s < t, and any w,w' such that | Re(rc) — a — [nt\c\ < e, 

| R e(w') — a — [nsjc| < e and \w — w'\ < 1. Then there exists a constant C = C(cr,T), depending 
only on a and T, such that 


E 


M [ntiAT( w ) - m \us\/\t{ w ') <C(a,T)c 4: \{t,w) - (s,w')\ 4 , 


where |(t,u>) — (s, ra / ) | denotes the Euclidean norm in 
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Proof. Clearly the l.h.s. of the above inequality is upper bounded by 


C 


E 


M \nt\AT^ W ) - M [nsj At(®< 


+ E 


M\ns\KT( w ) My nsjAT^T ) 


(19) 


for some absolute constant C > 0. We control the two terms separately. 

Set G 0 (w ) = G(w) — w, so that M/ c+ 1 (iu) — M^(w) = Go(Ffc(w) — i@k+ i) + c, and 

E(M k+ i(w) - M k (w)\cr(Qi, ... ,0*0) = ^ J G 0 {tk(w) - *i?)di? + c = 0 

for all k < [nt\ AT. It follows that {MkAT( w )) k<\ n t\ is a martingale. Moreover, we deduce from 
the estimates in ([3]), which now read 


|G 0 (in) + c| < 


Cc 


Re(w) — 5 ’ 


IG'oHI < 


Cc 


(Re(u>) — 5) 


2 ’ 


( 20 ) 


that |M ( fc+1 ) A7 ~(w) — Mk/\'r(w)\ = |G'o(ffc(u>) — i@k+ 1 ) + c| < Cc/a for all k < [nt\ . Using orthog¬ 
onality of the increments of (Mfc A 7 -), then, we find 


E 


m LT|at(t) - M [ ns JA t( W ) 


= E 


[nt J — 1 


[nt J — 1 

Y ( M (fe+1 )AtM - M kAr( w )) 

k=[ns J 


( 21 ) 


< 


C([nt\ - [ns\f Y E {\ M (k+i)AT( w ) ~ M kAT(w)f 
k= [nsj 


Cc 4 , 

< —o- t - s\ 

<7 8 


for some absolute constant C > 0 and n large enough. 


Let us now look at the second term in (19). For k < |_rasj A T set Mf. = M^{w) — Then 


{^kAT)k<\nt\ is a martingale, and by @ 


_ _ (J c 

\ M (k+i)AT ~ a 4at| = |G'o(ffe(u’ — i@k+ 1) - G' 0 (ffc(u’ / - *0fc+i)| < f J 2^\ w ~ w '\ + I^4atI) • 
This, together with orthogonality of the increments, yields 

It ®!- 1 B \ R®!- 1 

M \ns\AT 


E 


= E 


l/OJ x g\ " L 

Y {M(k+i)AT - MkAT) j < L ns J Y V(\M {k+l)AT -M kAT \ 

k =0 ' k =0 


< LnsJ 3 ^- Y (k-^'| 8 + IE(|MfcAr| 8 )) • 


[ns J — 1 


k =0 

It then follows from Gronwall’s inequality that 

14 Cc 8 


E 


Mi 


LnsJ AT 


< |_nsj‘ 


a 


16 


C[nsJ^c^ 

w — w'\ 8 e Tb < C(a,T)c 4: \w — w \ 


( 22 ) 


for n large enough and a constant C(a,T) depending only on cr and T. Putting (21) and (22) 
together, and using that \w — w /| 8 < \w — w’y since \w — w'\ < 1 by assumption, we find 


E 


M lnt\ATi w ) ~ M [ ns JA t( w ) 


+ E 


Mi 


[nsj AT 


/\|4 


< C\a, T)c 4 | (t, w) - (s,w')\ 


for n large enough and C'(a , T ) = max{C/cr 8 ; C(a, T)}, as claimed. 


□ 
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Kolmogorov’s continuity theorem now yields the existence of a random variable ^ > 0 such that 


sup 

\t~s\<ri 

w :| 'Re(w)—<r—[nt\c\<£ 
w':\ Ke(w') —a— [ns\c\<e 
\w—w'\<4:r) 


M\nt\KT( w ) - M \ns\AT{ w ') < M\(t,w) ~ (s,U>')| 1/16 , 


with E(..# 8 ) < C(a , T)c 4 , for C(a , T ) as in the statement of Lemma® Therefore we find 


sup 

|*—s|<T7 

w :| He(w)—cr— [ntjc|<e 
w':\ Ke(w')—a— [nsjc|<e 
\w—w'\<4:rj 


M lnt]Ari. w ) ~ M \ns\Ar{ w ') 


> vJc < 


> 


z/ 8 c 4 \ 




and sending first n —> oo and then r/ —> 0 we conclude that © holds. This proves tightness, and 
hence it concludes the proof of Theorem [4j 


5. The fluctuation process on J4? 

Notation. Recall that T = {z € C : \z\ = 1}, and set aT = {z G C : \z\ = a} for any a 6 M+. 
Let (C'(aT), || • Hoc) denote the space of continuous function on aT equipped with the supremum 
norm. Moreover, for a subset D of the complex plane, denote by the space of holomorphic 

functions on D. Whenever D = {|^| > 1}, denote Jff’(D) simply by Jf?. 

We have shown in the previous section that, for any fixed <r > 0, F° —> } :a as n —> oo in 
distribution with respect to the Skorokhod topology on the space D[ 0, oo) of cadlag functions from 
[0, oo) to C'(T). Note that, for any t > 0, the continuous function F%(t, ■) coincides with the 
restriction of the holomorphic function 

F n (t, z) = (log ^ ^ - [nt\ c) , (23) 

dehned for all \z\ > 1, to the circle e CT T. We show below that also the limit F a (t, •) can be in¬ 
terpreted as the restriction of a holomorphic random function F(t,-), dehned for all \z\ > 1, to 
e°T. Moreover, we provide an explicit construction of F, and prove that F n —> F in distribution 
as n — > oo, with respect to the Skorokhod topology on the space of cadlag functions from [0, oo) to 


Define, for all N > 1, a distance d/v on J4? by setting 

and let i/j) = 


d N ((f>,ip) = sup \(j>(z) - ip{z)\ A 1, 


N>1 


d N (<t>, VO 
2 n 


(24) 


Since (J4?, dm) is a complete separable metric space for all N > 1, this makes (J^, d) into a complete 
separable metric space (i.e. Polish space). It follows that, if Djp[ 0, oo) denotes the space of cadlag 
functions from [0, oo) to Jt? equipped with the Skorokhod metric <d#q then also (Djp[ 0, oo), d ,#f ) is 
a complete separable metric space. It is clear that each F n dehned in (23) is a random variable in 
Djfl 0,oo). 

Let us now describe the explicit construction of F. Let D = M/27rZ, and recall that if we set 
Sk(d) = e tk ^/\/2TT for fcgZ, then (ek)k&z forms an orthonormal basis (in short ONB) for L 2 (D) 
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with respect to the inner product (/, g) = Jff dd. On this basis, that we refer to as Fourier 

basis, the Poisson kernel Re (|^f) reads 


/ 1 I p^^ /T* \ _/ _ 

P 1 /r (tf) = Re ( 1 _ e ^ r ) = ^(eo + ^ 


r |fc| e fc 


fcez\{o} 

for any r > 1. Take two independent collections (fi k )kez and (P' k )kEZ of i.i.d. Brownian Motions on 
M, and denote by (A k ) kG i and (B k ) ke % the solutions to 


f d A k (t) = -\k\A k (t)dt + \/2d/3fc(t), 

|A*(0) = 0 


idB k (t) = -\k\B k (t)dt + V2d/3' k (t) 
(B k (0)=0. 


Then A k ,B k perform independent Ornstein-Uhlenbeck processes on M with invariant distribution 
jV(0, l/|fc|). Define formally 


WM)= Y, 

fee z\{o} 


A k {t) + iB k (t) 


(25) 


for (t, ?9) 6 [0,oo) x [—7r, 7r). Finally, for ( t,re m ) € [0,oo) x {|^| > 1} set 

1 r /U e 'M)/n 


Ht, re”) = -D (p 1/r » W(t, •)) (a) = 


\/27T 7-, 


Re 


Theorem 5 . F is a random variable in Cj?[ 0, oo) C Dj^[ 0, oo). Moreover, F n F as n ^ oo in 
distribution with respect to the Skorokhod metric d on -D^[0,oo). 


The question of how to make sense of the boundary values W defined formally in (25) is addressed 
in the next section, where we show that W can be rigorously defined as an Ornstein-Uhlenbeck 
process on a suitable infinite-dimensional Hilbert space. 

The rest of this section is devoted to the proof of Theorem [5} 

Lemma 9 . F is a random variable in Cjg>[ 0,oo). Moreover, F is Gaussian, and its restriction to 
[0, oo) x e°T agrees in distribution with F a defined in Theorem Q for any a > 0. 

Proof. The fact that F is Gaussian is true by construction. Fix any t > 0, and expand F(t, •) in 
Fourier basis, to get 


F(t,re ia )= ( e 0 + Y r ~ ]k]e ~ ika e k , Y 


k^O 


fc /0 


A k (t) + iB k (t)\ 


V2 


e k | = > r 

k^O 


y^ r -|fc| ( Ak(t) + iB k (t.) ^ ^ ika 


(d) 


r k [A k ( t ) cos ka — B k ( t ) sin ka] + i Y2 r k \Pk (t) cos ka + A k ( t ) sin ka] 


k> i 


k> 1 


(26) 


for any r > 1 and a £ [— vr,7r), where the last equality holds in law as stochastic processes, and it 
follows from the independence of the OU processes. This provides an almost surely convergent power 
series expansion for F(t,-) at all 2 with \z\ > 1, and hence it shows that F is a Gaussian stochastic 
processes taking values in Jf?. Recall from (24) the definition of the distances ( d]y),d on J#?. To 
see that F is continuous, we have to show that for all compacts of the form [0, T] it holds almost 
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surely that, if t n —> t E [0,T] as n —> oo, then d(fF(t n , •), J 7 ^, •)) -A 0, i.e. d/v(.T(f n , ')) ^ 0 

for all IV > 1. We have: 

d N (J r (t n ,-),J 7 (t,-)) = sup l-F^z)-.F(i,z)| = sup - .F(t,2r)| 

I^e 1 /^ |^| =e 1 / JV 

= sup I ^2 e~ k/N [(A k (t n ) + iB k (t n )) - ( A k (t ) + ?:B fc (f))]e* fca 

oe[-ir,7r) fe > 1 

< J] e " fc/JV (l^(*") ~Mt)\ + \B k (tn) - B k (t) l) . 

k> 1 


To show that the last term converges to 0 as n —> oo almost surely for all t n ,t E [0,T], then, it 
suffices to prove that the sequence of continuous functions gAt(t ) = YllLi e ~ k//N \Ak(t)\ converges 
uniformly on [0, T] as M —» oo. Indeed, we find 


OO OO 

sup \gM(t) - goo(t)\ = sup V] e" fc/JV |,4fc(f)| < V] e^ A:/iV sup \A k (t)\. 

i6[°,T] k=M+l k=M +1 iS[0,T] 


On the other hand Doob’s maximal inequality for the subnrartingale (e 2kt A k (t))t<T yields 


sup |A fc (t)| > e k/2N ) < P( sup {e~ 2kt A 2 k (t)} > e k/N ~ 2kT 

T 1 ! / V /r- fn T 1 


te[0,T] 


te[o,T] 


< 


E(^(T)) 




< 


1 


ke k/ N ’ 


so by Borel-Cantelli sup |klfc(t)| < e k ^ 2N for almost all k > 1, almost surely. This proves uniform 
te[o,T] 

convergence on compacts, and hence almost sure continuity of T. 


To conclude the proof, we show that T has the same covariance structure of J- a on the circle 
e a T, for arbitrary er > 0. Indeed, it follows from (26) that real and imaginary parts of J-(t, e CT+M ), 
a E [—7r,7r), are independent centred real Gaussian random variables, with 


E[(ReJ 7 (t, e CT+ * a )) 2 ] =E[(Im T(t,re a+ia )) 2 ] 


0 —2 ka 


(1 - e~ 2kt ) 


= log 


i _ e - 2 h+ CT ) 


k> 1 


1 — e 


— 2(7 


Moreover, expanding J r (s,e a ) as in (26), one further checks that for s < t 
Cov(ReJ 7 (t, e a+m ), ReJ 7 (s, e CT )) = Cov(lmJ'(t, e CT+M ), IraF(s, e a )) = Rc^log i _ e _ {t _ s) _ 2a+ia J , 
and 

/ 1 _ p-(i+s)-2cr+ia 

Cov(ReJ 7 (t, e CT+m ), Im.F(s, e a )) = -Cov(lmJ 7 (f, e^ 10 ), ReJ 7 (s, e a )) = Im( log 1 _ ^ {t _ s) _ 2a+ia 


By rotational invariance in the spatial coordinate, this is enough to conclude that J- and J- a have 
the same covariance structure, and hence the same law, on every circle of the form e°T, for arbitrary 
a > 0, as claimed. □ 


Proof of Theorem [5| For any N > 1, denote by Hjy the operator that maps a continuous function 
on e 1 / iV T to its holonrorphic extension to the outer region {|z| > e 1/W }. Moreover, let Z?jv[ 0, oo) 
and Djp N [ 0, oo) denote respectively the space of cadlag functions from [0, oo) to (C(e 1 / jV T), || ■ |joo), 
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and the space of cadlag functions from [ 0 , oo) to (J4?({\z\ > e 1//iV }), d]\r) , both equipped with the 
Skorokhod topology. Finally, let /n ■ Dn[ 0 ,oo) — > D^ N [0, oo) be defined for x £ Dat[0,oo) by 

f N (x) : [0,oo) -A- J^({\z\ > e l,N }) 
t H N (x(t,-)). 


We claim that the proof of Theorem [5] amounts to showing that the map /jy is continuous. Indeed, 
assume so. Then it follows from Theorem [4j together with the continuous mapping theorem (cf. [I], 
Theorem 2.7), that .fiviJ'n^) —> /n^ 1 ^) in distribution with respect to the Skorokhod metric on 
A^[0,oo). Moreover, recall from ( |26[ ) the definition of T. Then, since by Lemma [9] the stochastic 
process jri/ N a g rees i n i aw with T on [ 0 , oo) x e ly/Ar T, we deduce that the corresponding holomorphic 
extensions /n^ 1 ^) and T agree in law on [0,oo) x {|z| > e 1/W }. It follows that fN{^v/ N ) —> T 
in distribution with respect to the Skorokhod metric on Djg> N [ 0,oo). Since fN{^n N ) = T n on 
[0, oo) x {\z\ > e 1 /^}, with J~ n defined as in (23), this shows that T n -» T in distribution with 
respect to the Skorokhod metric on Djp N [ 0, oo). N being arbitrary, we conclude that T n — > T in dis¬ 
tribution with respect to the Skorokhod metric d ^ on D 0, oo), which is what we wanted to show. 


It remains to prove that for all iV > 1 the map /w is continuous from Thv[0,oo) to Djg> N [0, oo). 
Since fjy only acts on the spatial component, it suffices to show that the holomorphic extension 
map iL/v is continuous from (C'(e 1 / iV T), || • Hoc) to > e 1//iV }), c?tv) for all N > 1. Indeed, 

take ( (f> n )n, 4> hr C{e 1 ^ N T) and suppose that \\4> n — 4> ||oo —> 0 as n —> oo. Then we have: 

d N (H N ((j) n ),H N O)) = sup \(H N (j) n )(z) - (H N <f>)(z)\ A 1 = sup \(H N (j> n )(z) - (H N (j>)(z)\ A 1 

\z\>e 1 / N | z\=e*-f N 

= sup \cj>n(z) - 4>(z)\ A 1 < ||0„ - <j6||oo 0 

| z\=e 1 / N 

as 7i —y oo, where the second equality follows by applying the maximum principle, which is in force 
since (H^()) n ){z) -» 0 as \z\ —> oo by construction for all N, n > 1. This concludes the proof. □ 


6. The boundary process 


We have proved weak convergence to a limiting Gaussian process T taking values in the space 
of holomorphic functions on {|z| > 1}, of which we have provided an explicit construction. In this 
section we address the question of a rigorous definition of the boundary values W of J~, formally 


given by (25) 


6.1. Abstract Wiener Space construction. It is clear that we have no hope to define W(t, •) 


pointwise, since the formal series (25) diverges almost surely at each point. One could try, on the 


other hand, to make sense of it as a complex Gaussian process taking values in a suitable Hilbert 
space of functions on the unit circle. 

For convenience of the reader, we review below the construction of Gaussian random variables, 
and then Gaussian stochastic processes, on infinite dimensional Hilbert spaces. This will then be 
applied to rigorously define the boundary process W. Our presentation follows m Eg. 
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6.1.1. Gaussian random variables on a Hilbert space. 

Definition 6. An Abstract Wiener Space is a triple (H,B,p,), where: 

(i) (H,(- , ■)#) is a Hilbert space, 

(ii) (-B, || • ||#) is the Banach space completion of H with respect to the measurable norm || • ||# 
on H , equipped with the Borel a -algebra B induced by || • | #, and 

(iii) /x is the unique Borel probability measure on (B, B) such that, if B* denotes the dual space 
of B, then /a o </> _1 = J\f(0, ||^||#) for all (f> G B*, where (j) is the unique element of H such 
that 4>(h) = ((f>, h)n for all h G H. 

Note that (iii) also reads as follows: if X is a random variable on (B, B) distributed according 
to //, then <j>(X) ~ A/"(0, ||0||#) for all (j) G B*. In this case we sajQthat X is a standard Gaussian 
random variable on H. We refer the reader to HUES] for the definition of measurable norm on a 
Hilbert space, and for existence and uniqueness of such a measure fi. Here we will only need the 
following two properties. 

Facts. Let H be a Hilbert space with inner product (-, •)#, and let || • |||, = (•, •)#. Then: 

(a) any measurable norm on H is weaker than || • ||#, and 

(b) if T is a Hilbert-Schmidt operator on H. i.e. 

OO 

E llTe, HI, < oo , ( ei )? =1 ONB of (H, (-, -) H ) , 

i= 1 

then || T • || # is a measurable norm on H. 

It is worth to point out that, unless dim (H) < oo, the norm || ■ |||, = (•, ■)# is not measurable on 
H , so that || • ||# ^ || • ||#. When || • ||# is itself induced by an inner product, which will turn out to 
be the case in our construction, B is itself a Hilbert space. 

6.1.2. Brownian Motion on a Hilbert space. 

Definition 7. Let (H, B, /x) be an abstract Wiener space, and denote by /# the uniqucj^] probability 
measure on B , such that /# ° <p = A/"(0, t||^||#) for all 4> G B*. Finally, let C#[0,oo) the space of 
continuous functions from [0, oo) to B, equipped with the tr-algebra generated by the coordinate 
functions x i-A x(t). Then (cf. [5], pp.81-85) there exists a unique probability measure /x on C#[0, oo) 
such that, if IT is a random variable in C#[ 0 , oo) distributed according to /x, then the following 
hold: 

• IF(0) = 0 /x-a.s. 

• W has independent increments, 

• for any 0 < s < t, W(t) — W{s) is distributed according to /xt_ s . 

If a random variable IT on C#[0, oo) is distributed according to /x we call it (cylindrical) Brownian 
Motion on H. 

2 Note that, unless dimll < oo, a standard Gaussian random variable X on H does not take values in H, but only 
in the larger Banach space B. 

^Existence and uniqueness follow trivially by existence and uniqueness of /x. 
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Proposition 4 ([5], Proposition 4.3). Let (H,B,p) be an Abstract Wiener Space, and let (e k ) k be 
an ONB of H with respect to (•,•)#. If W is a Brownian Motion on (H,(-,-)h), then there exists 
a collection of i.i.d. real-valued Brownian Motions (/3k)k such that 

w(t) = ■£ f3k{t)ek, 

k 

where the above series converges in L 2 (B). 

Note that by setting t = 1 in the above result we deduce the following. 

Corollary 3. Let (H, B, p) be an Abstract Wiener Space. If X is a standard Gaussian random 
variable on H, i.e. X ~ p, then there exists a collection (Ak)k of i.i.d. jV(0,1) real random variables 
such that X = E Ak^k ■ 

k 

6.1.3. OU process on a Hilbert space. Having constructed a Brownian Motion W on H, one could 
then go ahead and define stochastic integration with respect to it. For a detailed account on 
the theory of stochastic integrals with respect to a Brownian motion taking values in an infinite¬ 
dimensional Hilbert space we refer the reader to [5], Chapter 4. Here we are only interested in a 
very special case, namely the one of deterministic integrands which diagonalise on the ONB (ek)k 
of H. Indeed, in this case the definition of stochastic integral with respect to W reduces to the one 
of a countable collection of stochastic integrals on M. 

Definition 8. Let ( H,B,p ) be an Abstract Wiener Space, and assume that the norm || ■ ||b on B 
is induced by an inner product (•, • )b , so that B is itself a Hilbert space. Let (ek)k denote an ONB 
for H, and let W(t) = Ylk Pk(t) e k be a Brownian Motion on H. Then, if 'P(f) = ^2 k a k (t)e k f° r 
some collection of continuous real-valued functions ( a k (t)) k such that Ylk fo° |nfc(t)| 2 di < °o, then 

/ T(s)dlF(s) := V f a k (s)d/3 k (s)e k , Vf > 0 . 

Jo k Jo 

6.2. The boundary process W. Let us now see how the above definitions read in the case of 
our interest. Recall that D = M/27tZ, and for /, g £ C°°(D ) with Fourier expansion /($) = 
Efcez A e fc(tf), = EfcezfffcOc(tf) introduce the inner product 

U,9)h ■= \ k \fk9k- 
k^O 

Let ~ be the equivalence relation on C°°(D ) which identifies two functions if they differ by a 
constant, and set H S (D) := C°°(D)/ ~. We identify each equivalence class of H S (D ) with its 
representative having zero average on D. Denote by H the Hilbert space completion of H S (D ) with 
respect to (•,•)#. We seek to define an abstract Wiener space ( H,B,p ) for a suitable choice of 
measurable norm on H. 

By property (b), in order to construct a measurable norm || ■ ||# on H it suffices to find a Hilbert- 
Schmidt operator T on H, and set || • ||b = ||T- \\h- In fact, we are going to construct a one-parameter 
family of such operators, and hence of measurable norms on H . 
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For any a G K, let (—A) a be the operator that acts on L 2 (D) functions by multiplying the Fourier 
coefficients by \k\ 2a , that is 

(—A) a ( Y f k e k ) (0) := Y \k\ 2 a fke k (V). 

' fcez ' o 

Set 

U a = {/ G L 2 (D) : (-A)“/ G L 2 ^)}/ ~, 

and equip with the inner product ( f,g) a ■= ((—A) a /, (— A) a g). Then (H a ,(-,-) a ) is a Hilbert 
space, and in fact it is the Hilbert space completion of H S (D ) with respect to the inner product 
(-,-) a . Denote by || ■ || a the norm induced by (-,-) a , and note that if a > b then || ■ ||& < || • || 0 , 
so that H a C Hb- It follows that, whenever b < a, || ■ ||& is well defined on H a , and moreover 
\\f\\b = ||(—A) 6_a /|| a . In order for || • ||j, to also be measurable on H a , then, it suffices to show that 
(—A) b ~ a is a Hilbert-Schmidt operator on H a - To this end, note that if f k = (—A)~ a e k for k / 0, 
then (f k ) kj i o is an ONB for Ha- Moreover, 

£ ll(-A)‘-“AllS = £ ((-A)Vfc,(-A)‘A) = £ ((-A)‘-“e t , (-A)*-%) =2£**<"-“>, 

k +0 k^ 0 k +0 fc>l 

which converges if and only if b < a — 1/4. We have therefore proved the following. 

Proposition 5. Whenever b < a — 1/4 the norm || • ||& is measurable on 

Note that ( f,g)n = ((—A) 1 / 4 /, (—A) 1 / 4 #) by definition, so (—A) -1 / 4 provides a Hilbert space 
isomorphism between the spaces H and Hi/ 4 , that when convenient we identify. It follows that in 
order to get a measurable norm on H it suffices to have a measurable norm on Hi/ 4 . By Proposition 
[5j any norm of the form || ■ ||_ e for £ > 0 will do. 

Let B a denote the Borel cr-algebra on H a for the norm |j • || a . Then for any e > 0 there exists 
a unique Borel probability measure /i_ £ on (H- e ,B~ e ) such that (iii) in Definition [6] holds, i.e. 
Us o (p - 1 = jV(0, \\4>\\h) for all 4> continuous linear functionals on H- e , where 0 is the unique 
element of H such that (j>(h) = ((f), h)u for all h G H. It follows that (H, H- e , H- e ) is an Abstract 
Wiener Space for any e > 0. 


Definition 9. A random variable X is said to be a standard Gaussian on the Hilbert space H = 
H- 1/4 if X ~ as random variable on H- e , for all e > 0. 


Remark. This definition is consistent, meaning that if —£1 < —£ 2 , so that H- e 2 C H - ei , then by 
property (iii) of Definition [6] the restriction of /j_ £1 to H - £2 coincides with ^_ £2 . 


To conclude this abstract construction, we point out that X as above has a very simple expansion 
in Fourier basis. Indeed, if f k = (—A) 1 / 4 efc = e k />/\k\, then ( f k ) k is an ONB for H- 1/4 — H. It 
then follows from Corollary [ 3 ] that there exists a collection {A k ) k& i i.i.d. real AA(0,1) random 
variables such that 

X(tf) = Y, AMO) = Y -4=f M) ■ 

k^O k^O v l fc l 

Moreover, in light of Proposition [4] it is now easy to construct a Brownian motion W on H: 
simply take a collection of i.i.d. real-valued Brownian Motions (Pk)kez, and set 

WM) = Yj k (t)f k (0) = Y , 


k^O 


k jL 0 VW\ 
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where the above series converge in L 2 (%_ £ ), for any e > 0. From W we construct an Ornstein- 
Uhlenbeck process on H as follows. Let 'L : [0, oo) —>• B be defined by 


*(i) = E VW\e- ]klt fk = E e ~ We * > 

fc/o 

and set W(t) := /q — s)dW(s) for all f > 0. Then by Definition [i] we have 


Vi>M) = E[ f e-^- s) dp k (s) 

k +0 


efc(tf), 


so that the coefficients of W perform i.i.d. OU processes on M. Take an independent copy W of 
W. Then we have 


W(t,-) + iW(t,-) 

V2 


sr' f A k (t) + iB k (t)\ 

L(—^ 

b-£C\ V ^ 


for i.i.d. OU processes (. A k ) k , {B k ) k on M, which equals the r.h.s. of (25). This provides a rigorous 
construction of the boundary process W of the limiting holomorphic field T as an OU process on 
Ti-e, for any e > 0. 


Remark. Note that, by linearity, 


dW(i, •) = - 


Ei fc i 

k^ 0 


Ak(t) + iB k (t)\ 


V2 


)e k (-) 


dt + V2 'y ^ 

k^O 


d Pk(t) + id/3' k (t) 

V2 


efc(-) 


where we have denoted by (/3k), (/3' k ) the i.i.d. real Brownian Motions driving the OU processes 
(Ak), (Bk). This shows that W solves the Fractional Stochastic Heat Equation 


dW(t, •) = -(-A) 1/2 W(f, •) + y/2 4£(t, •) 


for £ complex space-time white noise on the unit circle T. 

Remark. We point out that the convergence result J- n — > T of Theorem [5] can be interpreted 
as a convergence result for the corresponding boundary values, seen as distributions acting on a 
suitable space of test functions. More precisely, let W n denote the boundary values of J- n , so that 
W n (t,!?) = /F n (t,e™) for (i,i9) 6 [0, 00 ) x [— tt,tt). As space of test functions we take 

{<p £ C 00 ^) : ip = P r * ip for some i/j € C 00 ^) and some r > l} . 

For each such ip, we have: 

(W n (t,-),<p) = (' W n (t,-),P r *4>) = (P r *W n (t,-),ij>) = (/F n (t, re 1 '), ip) (B(t, re 1 '), ip) 


as n —> 00 in distribution, as continuous functions on T. It would be interesting to understand if 
such convergence holds for a larger class of test functions, and ultimately as stochastic processes 
taking values in T-i- e for e > 0. 
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7. Further developments 

As mentioned in the introduction, arguably the main open problem in this area is to obtain 
rigorous results on (the non-regularised version of) HL(a) models with a > 0, with a particular 
interest for a E [1,2]. In this last section, on the other hand, we would like to collect some open 
questions for the a = 0 case, which has proven to be already very interesting from a mathematical 
point of view, and it has the advantage of being more tractable than the case a > 0 due to its 
intrinsic i.i.d. structure. 

Maxima of the cluster boundary. Once the question of fluctuations of cluster boundary has 
been settled, one is led to wonder about the asymptotic behaviour of the maxima of the cluster. 
For n > 0, set m n = sup{|z| : z E K n }. It follows from the analysis carried out in [25] that 
| m n — e cn \ -A 0 almost surely as n —loo, c A 0 and nc — > t for some t > 0. What about fluctuations 
around this deterministic behaviour? 

Connections with IDLA. Let us recall the definition of a different, discrete model of aggregation 
on the plane. Consider the lattice Z 2 , set Kq = {(0,0)} to be the initial cluster, and define 
K\ c K ‘2 C A 3 ... recursively as follows. At each step n > 1, start a random walk X n = (X n (k))k£^ 
(independent of everything else) from the origin, and, if T n = inf{7c > 0 : X n (k) ^ AT n _ i} denotes 
the first time X n exits the cluster K n _ i, set K n = K n _iU{X n (T n )}. This grows an increasing family 
of connected clusters (K n ) n > o C Z 2 , whose dynamics is usually referred to as Internal Diffusion 
Limited Aggregation (in short IDLA) on Z 2 . IDLA was originally introduced by Meakin and Deutch 
in [23] . and its large n behaviour is by now well understood in any dimension. We focus here on 
the 2-dimensional case, which we have described above. It is shown in [20] that K n ~ B r ^{ 0) 
for n large enough, where B r ( 0) denotes the Euclidean ball of radius r centred at the origin, and 
r(n) = sjn/ir , so that 7 rr 2 (n) = n. In mm fluctuations around this asymptotic spherical shape 
are analysed. In particular, the authors show that the maximum fluctuations of K n are O(logn) 
for large n. The average fluctuations are then studied in [16] . and are shown to be given by the 
restriction of a variant of the 2-dimensional Gaussian Free Field, so called Augmented GFF , to the 
unit circle. If, on the other hand, one considers IDLA on the 2-dimensional cylinder (Z/IVZ) x Z + 
rather than the whole plane, then (cf. |17| 1 the fluctuation field coincides exactly with the restriction 
of the whole plane GFF to the unit circle, i.e. our limiting fluctuation field Woo- One is then led 
to wonder whether the connection between HL(0) and IDLA clusters goes beyond having the same 
scaling limit and fluctuations. 

Log-correlation and branching structure. Log-correlated Gaussian fields appear to arise in 
correspondence with underlying branching structures. To make this statement more precise, let us 
give some examples. 

Branching Random Walk. For N E N+ consider a binary tree Gn = (Vn,En) of depth N, and 
assign to each edge e E En a standard Gaussian random variable X e in an independent fashion. 
Let Ln denote the set of leaves of the binary tree (so that |Ljv| = 2^), and for each v E Ln denote 
by Jv the unique path from v to the root. Then set 

Yv=Y 1 X e ■ 

eS7„ 
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It is easy to see that, when the set of leaves is seen as embedded in the unit interval [0,1] (i.e. 
the leaves are identified with 2 N equispaced points in [ 0 , 1 ]), the correlation between two leaves 
v,v' E Ln decays logarithmically with their Euclidean distance \v — v'\, for all N. In other words, 
at each level N > 1 the random field ( Y v ) v ^l n is a log-correlated (Gaussian) field. 

Branching in IDLA. A branching structure is also present in the IDLA model described above. 
Indeed, one could imagine to colour, each time a new particle is added to the cluster, the last edge 
along which the correspondent random walk has jumped in order to exit the cluster. More precisely, 
define the set of coloured edges recursively as follows. Given the cluster /\ n _i with n — 1 particle, 
start an independent random walk X n from the origin and, if T n is the first exit time of X n from 
K n -i, colour the edge ( X n (T n — 1), X n (T n )). It is easy to see that the set of coloured edges forms a 
tree. Moreover, if one considers two points on the cluster boundary (which is asymptotic circular), 
then it follows from the work in m that the correlation between fluctuations from circularity at 
the two points decays logarithmically with their distance. 

Branching in HL growth. The fact that a branching mechanism arises in the growth of HL(0) clus¬ 
ters was discovered by J. Norris and A. Turner in [25]. Think for simplicity of the particle P as 
being a slit, i.e. P = [1,1 + 5] C C. Then each time a new particle arrives it gets attached to 
exactly one particle in the cluster, at exactly one point, almost surely. Call this parent particle. 
This identifies an ancestral structure in the cluster. Given any finite set of points, it is shown in 
[25 ] that the corresponding ancestral lines converge to backwards coalescing Brownian Motions, 
i.e. Brownian Motions which evolve independently backwards in time until they coalesce, at which 
point they merge. This result can be rephrased by saying that the cluster boundary converges, in a 
finite-dimensional sense, to the so called Brownian Web. Here the underlying (forward) branching 
structure is then really a backwards coalescing structure of the cluster boundar}j^J but nonetheless 
one obtains log-correlated fluctuations. 

In all the above examples log-correlated Gaussian fields arise as fluctuation fields of several 
different models, all featuring some kind of underlying branching structure. We believe it would be 
very interesting to understand to what extent this is a general phenomenon, and, if so, how robust 
it is with respect to variations of the branching mechanism (e.g. introduction of correlation between 
branches). 


Appendix A. Local fluctuations 
Our main result on local fluctuations is the following. 

Theorem 10. Pick any t > 0, and let z = e m+a , w = e a for some a E [— 7 r, tt), a > 0. Define 
e = (5 2 / 3 log(l/5) as in Theorem [jj and assume that a — > 0 as c —> 0, with a 3> yfe. Then, as 
n oo, nc t and a —> 0, it holds: 

+ log+l-ncj 

\ v cl °g(^) v cl °g(^) / 

3 Note that the dynamics of backwards coalescing Brownian Motions is different from the one of forward branching 
Brownian Motions, since for example in the latter case paths can intersect after branching. 









FLUCTUATION RESULTS FOR HASTINGS-LEVITOV PLANAR GROWTH 


29 


in distribution, where (JVi, A/2) is a random vector with centred complex Gaussian entries, and 
covariance structure given by 

/ 1 a \ 

~ 1 PT _ 

1+Q! 2 / 


E(MM) = ^ 

\ l+o; 2 


for a = lim CT _>.o ^ G [0, 00 ] (with the convention that 1+ 1 q2 = = 0 when a = 00 ). 


Theorem 


10 


place of 


follows by the same arguments that lead to Theorem 


considering now , ' m 


\/ lo S^) 


Under the assumption a S> \fe Lemmas |3|4| still apply, so that Theorem |2j holds. 
Together with Proposition [2j this shows that in the limit as cr —> 0, n —> 00 and nc —> t, real and 

log ( /TT ) ~ nc 

imaginary part of- , em+<T — are asymptotically i.i.d. centred Gaussians, with limiting variance 


given by 


i/ cl °g(^) 


lim 


1 


ff ^ 0l °g(^) L 2?r 


ra+t 


Re 


3 — ifl+x 


+ 1 


g— ifl+x _ ^ 


dddx — t 


= lim- . -1 . 

log(^) 


log 


1 _ g-2(cr+i) 


1 — e 


—2a 


= 1 . 


For two point correlation we reason as in Section 4.1 to gather that the limiting covariance between 
. . 

and Ref 


Re 


log ■ 


e ia+cr 


lim 


\/ c 1 °s( 2 f) 

1 




V clo g(2F) 


is given by 


1 


r<r-\-t 


a ^° !og(^) L 27r 

= lim 


' — i($— a)+x 1 I 

Re | r- ) Re 


cr—>-0 


lo g (£) 


log 


g— i(P—a)+x J 

1 _ g—2 (a+t)+ia 


e -ii)+x + X \ 


dildx — t 


^ g—2 cr+ia 


= lim 

cr—>-0 


g— ifl+x _ ]/ 

log(l + e 4cr - 2e~ 2a 


cos a) 


2 log 2 a 


1 + a 2 


whenever a/2a -A a G [0, 00 ], where the last equality is obtained by Taylor expanding around 
a = 0, a = 0. The same holds for imaginary parts correlation. 


Finally, the asymptotic covariance between Re 


log 


$n( 




by 


lim 


1 


1 


rcr-\-t 


lo g(^) Y 1 ^ J a 

= lim 


e —i(jd—a)+x _|_ 1 
, ^ ( e -i(#-a)+x _ l 


Im 


cr—>-0 


lo g(£) 


Arg 


\ _ g—2 (a+t.)+ia 
2 g— 2cr+m 


V cl °g( 2F) 

e -iP +x + 1 
g —ii)+x _ 

= lim 


^~ nC ) and lm( los pt±? c ) 


is given 


o-->0 log 2(7 


diMx — t 


arctan 


sin a 


cos a — e 


2 a 


a 


1 + or 


whenever a/2a —> a G [0, 00 ], where the last equality is obtained by Taylor expanding around 
a = 0, a = 0. This concludes the proof of Theorem |10[ 


Appendix B. Proof of Lemma \7\ 

Fix any k < |_nsj and set for simplicity Z\ = Z k ^ nt j (a), Z2 = Z k ^ ns j (0), W\ = e ia + CT +(Ynti-k)c^ 
W 2 = e a +(Y ns i~ k ) c _ Moreover, introduce the functions g^(z) = Re( log F ^L i of' > ), h$(z) = cR.e( ^_‘^j 
Then E(X fc) |_ nt j (a)X k ^ ns ) (0)| = 2 k T 9 ^Z 1 )g^(Z 2 )dd - c, and we have to show that 


T/ OT (z l)w (z 2 )dtf-T 


^(bFi)^(VF 2 )di? 


< 


C(t)ce 


(cr + (|_nsj — &)c ) 3 ’ 


(27) 
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Trivially, the l.h.s. is bounded above by 

2 ^: f (M^)| ■ \g^Z{) -M^)l + IMWOI • |^(Z 2 ) -M^2)|)di?. 

We bound each term separately. Recall the definition of E(rn,e), from which it follows that 
max {| Z Untl (0)| V |Z fc>M (t?)|| < e "+(L"tJ-*)*(! + 2e) < 2e t+ 2 = C(t ), 

min {|Z* iLntJ (i?)| A |Z fc)LrwJ (0)| - l} > e -+(M-^(i _ 2e) > a + ([ns\ - k)c 

tfe[-7r,7r) 1 L J J 2 

as long as n is large enough and a /§> e. We combine the above estimates with the bounds in 
Corollary [2j to get 

Z 2 + 1 


\g&{z 2 )\ < 


, F(e~^Z 2 ) e^Z 2 +1 

log- 7 ^-C- 


< 


e-«Z 2 

C(t) C 3 / 2 


5 “*% - 1 


+ C 


3 “*% - 1 


— (|z 2 | — i) 3 |Z 2 |-1 


+ 


C(t)c 


< 


2C(t)c 


(28) 


(a + (|_nsj — fc)c) 3 cr + (|_nsj — &)c cr + ([nsj — k)c 
for n large enough. Similarly, we find 


IMWi)| <c i + 


|Wi|-l 


<c 1 + 


< 


C(t)c 


<7 + (|_nsj — fc)c / cr + (|_nsj — k)c 


(29) 


on the event E(m, e) for, say, C(t) = t + 3 and n large enough. 

In order to bound the remaining terms, observe that, by definition, on E(m,s) we have 

max { \Z K[ns ^) - e ^+(Ln S j-fc)c| y | Z (l?) _ e W+«T+(L«tj-fc)c|| < c- (i)e> 

17G[—7T,7r) ^ J 

from which we get |Zi — Wi| < C(t)s and \Z 2 — W 2 \ < C(t)e. Combining this with Corollary [2] we 
finally obtain 


|<?(^i) - h(Wi)| < | 5 (Zi) - h(Z 1 )| + |h(Zi) - /r(Wi)| < 


Cc 3 / 2 |Zi| 2 


+ 


2c|Zi — W\\ 


{a + {[ns} _ k ) c y i\Zi\ — 1 )(|Wi| — i) 


< 


C(t)c 


,3/2 


+ 


C(t)c£ 


< 


2 C(t)ce 


(cr + (|_nsj — /c)c) 3 (<t + (LnsJ — fc)c) 2 (cr + (|_nsj — k)c) 2 

for n large enough. Similarly one shows that the same bound holds for \g{Z 2 ) — h(W 2 )\. Putting 
now this together with (28) and (29) gives (27). 

The second statement of the lemma follows by the same arguments, and the proof is omitted. 
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